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We study the spin dynamics of carriers due to the Rashba interaction in semiconductor quantum 
disks and wells after excitation with light with orbital angular momentum. We find that although 
twisted light transfers orbital angular momentum to the excited carriers and the Rashba interaction 
conserves their total angular momentum, the resulting electronic spin dynamics is essentially the 
same for excitation with light with orbital angular momentum I — +|Z| and I = — |Z|. The differences 
between cases with different values of |Z| are due to the excitation of states with slightly different 
energies and not to the different angular momenta per se, and vanish for samples with large radii 
where a fc-space quasi-continuum limit can be established. These findings apply not only to the 
Rashba interaction but also to all other envelope-function approximation spin-orbit Hamiltonians 
like the Dresselhaus coupling. 

PACS numbers: 78.20.Bh, 78.20.Ls, 78.40.Fy, 42.50.Tx 


Light with orbital angular momentum (0AM), referred 
to as twisted light, is a relatively new field of research 
which has become increasingly popular [1-17] since Allen 
et al. showed how twisted light beams can be easily gen¬ 
erated from conventional laser beams [18]. Recently, the 
theoretical foundation of the optical excitation of solids 
and nanostructures with twisted light has been estab¬ 
lished [19-27], and experimental studies with twisted 
light on semiconductors have been carried out [28, 29]. 

One motivation for such studies is the prospect of us¬ 
ing the large amounts of angular momentum that twisted 
light can carry in order to control the spin dynamics of 
electrons, thus adding a flexible tool to the active field of 
spin control [30-39]. In this context two different mech¬ 
anisms need to be distinguished. First, angular momen¬ 
tum as well as energy selection rules can lead to selective 
optical excitation of carriers with a preferred spin direc¬ 
tion. This mechanism enables fast spin-selective prepara¬ 
tion of states during the photoexcitation process and has 
recently been studied for strongly confined systems such 
as quantum dots [27] and quantum rings [23]. Secondly, 
the spin-orbit interaction—like the Rashba [40] and Dres¬ 
selhaus [41] couplings in semiconductor structures—is ex¬ 
pected to couple the 0AM of carriers transferred from 
the twisted light [19, 22] to their spin degree of freedom. 
This would provide a slower carrier spin control which 
would be dynamical and would remain active after the 
twisted light pulse. 

In this Letter, we study the spin dynamics of carriers 
in semiconductor quantum disks and wells excited with 
twisted light taking into account the Rashba spin-orbit 
interaction. Our central finding is that, rather unexpect¬ 
edly, the spin dynamics of the photo-excited electrons 
differs only slightly after excitation with light with and 
without 0AM in the limit of large quantnm disks, be¬ 
coming insensitive to the 0AM content of the twisted 


light beam for extended quantum wells. This result is 
consistent with the outcome of recent experiments which 
did not show traces of the 0AM transferred from twisted 
light to bulk GaAs in spin-resolved photoemission mea¬ 
surements [29]. 

Analytically, we find that the Rashba interaction, while 
conserving the total angular momentum of the electrons, 
has matrix elements which are independent of the 0AM 
quantum number in the (k-space) quasi-continuum limit. 
As a consequence, the induced spin dynamics is almost 
identical, in particular, for twisted light with components 
of the 0AM in the growth direction I = -|-|Z| and I = —\l\. 
This finding can be generalized to all possible effective 
spin-orbit interactions stemming from a lattice-periodic 
potential in the envelope-function approximation, e. g. 
the Dresselhaus coupling. From this we conclude that 
the dynamical spin control mechanism as analyzed here 
can only be effective for small quantum disks and other 
strongly confined systems. 

The discussion of the optical excitation of electrons 
with twisted light is especially clear when a basis of 
cylindrical states is chosen [22]. The wave functions of 
these basis states are expressed in cylindrical coordinates 
{r, (j), z} as 

tpbmi^ir, (t>, z) = Afmi^Jm{kTni^r)e'-^'^^biz), (1) 

where Jm is the m-th Bessel function, $t,(z) is 
the z-envelope of the b subband (the band index 
includes the spin quantum number) and Afmu = 
is the normalization factor. For a 
circular quantum disk with radius R, width L and growth 
direction z, the boundary conditions ipbmi^{R,4’j z) = 0 
are satisfied, if kmi' = Ura,vlR where Ura,v is the ly-th 
zero of the m-th Bessel function. Note that ipbrniy is an 
eigenstate of the z-component of the envelope 0AM oper¬ 
ator with eigenvalue hm and kmi/ determines the kinetic 
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FIG. 1. Diagonal elements of the density matrix pcmvcmv after an excitation of a cylindrical quantum disk with radius R with 
a box-shaped pulse of length t = 5 ps and orbital angular momentum I of the light, (a) and (b) show the occupations of a 
R = 1.2 pm quantum disk with I — 0 and I = 5, respectively, (c) and (d) display the occupations for a larger i? = 5 pm disk 
with I = 0 and I = 5. For a better comparison, the matrix elements are plotted against the energies tcm/i instead of the indices 
p and the absolute values of the occupations are rescaled. 


energy of the state ipbmui since in a parabolic band b 
with effective mass the energy of the state is given 
by Chmu = Note also how the precise lo¬ 

cation of the energy eigenvalues is given by the zeros 
of the Bessel functions; this detailed information will be 
“smeared out” in the limit i? —oo as the allowed values 
of k become a quasi-continuum. For the sake of simplic¬ 
ity, we restrict our discussion to a case where only spin 
degenerate conduction (6 = c) and heavy-hole (h = v) 
bands are considered. 

The matrix elements of the twisted-light-matter inter¬ 
action Hamiltonian Hj (in the dipole approximation for 
only the z-component) in the cylindrical basis states was 
derived in Ref. 22: 

{cUl V \}r [= ^cvi/' um' ^ — h 

g 

rUe 

R 

X J *dv JjYi'{kjfif^ (2t)) 
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where e and me are the electron charge and (bare) mass, 
Aq is the field strength, uj is the light frequency, q\\ the 
in-plane and the growth-direction part of the light 
wave vector, is the electron wave vector in the growth 
direction, I the 0AM of the light, is the light polariza¬ 
tion vector and Pcv is the dipole matrix element between 
heavy-hole and conduction band states. Note that 
contains spin selection rules. Let us consider the case 
where due to excitation with circularly polarized light 
only spin-up electrons are excited. 


In Ref. 22 equations of motion were presented for 
the density matrix under the influence of twisted light 
switched on with constant amplitude at t = 0. In the 
low-excitation limit, i.e. initially empty conduction and 
filled valence bands excited with a moderate light field 
so that the occupations can be well approximated by an 
expansion up to second order in the field strength, we 
find from Eqs. (16) and (17) of Ref. 22: 
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Thus, the optical excitation yields only diagonal ele- ments of Pckmc'k'm' with respect to m, i.e. only states 
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with a defined envelope 0AM are excited. Also, we find 
from Eqs. (2) that for every electron with 0AM m a hole 
with 0AM m — I is excited. From this we can conclude 
that the total envelope 0AM induced in valence and 
conduction band together is = hlNf, where Ne is the 
number of excited electrons or, equivalently, holes. For 
very short pulse times t, the diagonal elements of the 
density matrix from Eq. (3) are given by 

~ ^ ^ 1^^ (—I 


Note that in Eq. (4) no information about the band struc¬ 
ture is contained. In particular, the formula for the va¬ 
lence band occupations is the same as for the conduction 
band. Thus, in this limit, the total angular momentum 
is distributed symmetrically to heavy-hole and conduc¬ 
tion band, i.e. the total envelope 0AM in the conduction 
band is = H^Ne- However, for longer pulse durations, 
the energy selection becomes important leading to an in 
general different value of 

Figure 1 shows the diagonal elements of the density 
matrix Pcmvcmv after an excitation with pulse duration 
t = 5 ps with circularly polarized light with 0AM 1. 
The central frequency of the light pulse was chosen to 
be resonant with the band gap. The effective masses 
were m* = 0.067 me and = 0.45 me for conduction 
and heavy-hole band, respectively. The oscillatory struc¬ 
ture of the occupations along the energy axis can be at¬ 
tributed to the finite pulse duration via the energy-time 
uncertainty relation. Along the m axis, there are also 
oscillations in the occupations. Since their frequency de¬ 
pends strongly on the radius R of the sample and they 
get smeared out for large values of R, we attribute these 
oscillations to finite size effects. Note that in Fig. 1(b), 
where the occupation for light with Z = 5 is plotted, the 
states with the 5 lowest values of m for every energy shell 
are empty (seen more clearly at low energies), since there 
are no valence band states which satisfy the condition 
m' = m — I oi the matrix element in Eq. (2a). Figures 
1 (c) and (d) show that for the larger R = 5 /rm quantum 
disk, the difference between the occupations after I = 0 
and I = 5 excitations diminishes visibly. 

Now, we focus on the spin dynamics after the optical 
excitation. We study the effects of spin-orbit coupling 
mechanisms, considering for concreteness the Rashba 
Hamiltonian [40], Hu, which is usually the dominant 
mechanism in quasi-two-dimensional systems. In order 
to better work with the cylindrical states given in Eq. (1), 
we switch from the usual cartesian-coordinate expression 
of Hji to its expression in polar coordinates: 


Hr =haR{kyax - kxCTy) = HaH{s^d -s 9+), (5a) 
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FIG. 2. Spin dynamics after excitation with twisted light 
with orbital angular momentum I = {—5,0, 5} and quantum- 
disk radius R = 1.2 pm and with I — {0,5} and R = 
{5 ^m, 10 /im}. 


where aR is the Rashba coefficient, ai and are the 
Pauli matrices and spin raising and lowering operators, 
respectively. With the relation 
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it is straightforward to calculate the matrix elements of 
Hr with respect to the cylindrical states: 


O P Ic / I 

{cmv\d \cmv) = -, (7a) 


(c'm'u'lHjilcmu) = 
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It can be seen from the form of the Rashba Hamiltonian 
in cylindrical coordinates that an electron with spin-up 
(down) and envelope 0AM m flips to a state with spin- 
down (up) and 0AM m -I- 1 (m — 1). In this sense, d^ 
can be regarded as raising and lowering operators in m. 
If Hr is applied a second time, the electronic state is 
transferred back to the initial spin and 0AM state, while 
a change in is possible. Note that the sum J = m -I- s 
of the envelope 0AM m and the spin s is conserved by 
the Rashba Hamiltonian. 

Having derived the matrix elements of Hr in cylin¬ 
drical coordinates, it is straightforward to calculate nu¬ 
merically the time evolution of the density matrix, where 
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the initial conditions correspond to the final occupations 
generated by optical excitation with light with 0AM /, 
illustrated in Fig. 1. The resulting dynamics for the 
total conduction-band spin is shown in Fig. 2. We 
show results for three different values of the disk radius, 
R = {1.2, 5,10}/im. As in the case of optical excita¬ 
tion with light with zero 0AM, the Rashba interaction 
leads to a dephasing of the initial electron spins. Since 
for small disks only a finite number of states contributes 
noticeably to the dynamics, oscillations are found which 
do not cancel completely so that for long times the to¬ 
tal spin reaches a non-zero value. Note that the curves 
for excitation with I = 5 and / = — 5 coincide (shown 
for R = 1.2/im). This unexpected result shows clearly 
the insensitivity of the spin dynamics, in the presence of 
the Rashba spin-orbit interaction, to the content of 0AM 
transferred from twisted light to the electron gas. For the 
same quantum disk, an excitation with I = 0 produces 
spin dynamics slightly different from the I = ±5 excita¬ 
tion, but this difference decreases for larger radii. This 
tendency can be seen by comparing the excitation with 
I = 5 and I = 0 for the three different values of R used 
in Fig. 2. 

To understand this finding, it is useful to analyze the 
case of an infinitely extended quantum well, obtained 
letting i? —)■ oo. In this limit, the discrete kmi/ become 
continuous and the eigenstates can be written as 

'ipbkm{r,(j>,z) := (8) 

with the orthogonality relation {bkm\b'k'm') = 
Sbb'Smm'S{k — k'). Using that Bessel functions satisfy 

OO 1 

Jdr rJmikr)Jm{k'r) = j6{k - k'), (9) 

0 k 

the corresponding matrix elements become 

{c' k'm! \d^\ck'm) = z^kSm' ,m±i5{k — k') (10a) 

{c'k'm'\Hii\ckm) = hank 6{k — k')x 

). (10b) 

Thus, in the quasi-continuum limit, the prefactor of the 
Rashba-interaction depends only on the energy of the 
state via k but not on m. The spin dynamics is there¬ 
fore a precession of the electron spin with a fc-dependent 
frequency which, for a given k, is the same for all differ¬ 
ent values of m. Since the effect of the excitation with 
twisted light was mainly that states with different m are 
excited, it is now easy to see why, for extended systems, 
the spin dynamics due to the Rashba Hamiltonian is al¬ 
most the same for excitations with light with and without 
0AM. 

It is noteworthy that this statement is also true for all 
effective Hamiltonians with a microscopic origin in the 
lattice-periodic crystal potential such as the Dresselhaus 


[41] spin-orbit coupling when treated in the envelope- 
function approximation: For a lattice-periodic potential, 
the solutions of the corresponding Schrddinger equation 
are given by the Bloch theorem as ■i/'(r) oc e®'^''u„k(r) 
with lattice-periodic Bloch function Mrik(r), band index 
n and wave-vector k. The envelope-function approxima¬ 
tion consists of integrating over the plane-wave part of 
the wave function yielding an effective Hamiltonian [42] 
Hgff for the M„k which is diagonal in k and the matrix 
elements can be written as a power series in k. The 
resulting effective Hamiltonian can be rewritten by de¬ 
composing kx and ky in terms of d~^ and d~, as done 
in Eqs. (5). Thus, the dependence of the matrix ele¬ 
ments of Hgfi in cylindrical states on m is of the same 
character as for the Rashba Hamiltonian and vanishes in 
the quasi-continuum limit. For systems with finite size, 
however, a weak dependence on m can be found due to 
the TO-dependence of the possible fc-values in the prefac¬ 
tor of the Rashba Hamiltonian in Eq. (7b). This means 
that, e.g. for small quantum disks, where the energy sep¬ 
aration between the discrete cylindrical states becomes 
important, the 0AM of the exciting light can influence 
the spin dynamics significantly. 

In conclusion, we have shown that, although the orbital 
momentum of light can be transferred into the envelope 
orbital angular momentum of electrons, the usual solid 
state spin-orbit interactions, such as Rashba and Dres¬ 
selhaus interactions, do not couple the envelope orbital 
momentum of the carriers to the spin degree of freedom 
in such a way that a significant difference in the spin 
dynamics after excitation with light with and without 
orbital momentum is found in large extended systems. 
This finding can explain that in recent experiments [29] 
no influence of the light orbital angular momentum on 
the spin polarization was found. However, for cylindri¬ 
cal quantum disks with small radii, the discreteness of 
the states plays an important role so that the spin dy¬ 
namics indeed depends on the orbital momentum of the 
light. Nevertheless, also for small systems, the spin dy¬ 
namics after excitation with orbital momentum I = \l\ 
and I = —1/| are very similar, in contrast to optical exci¬ 
tation with opposite circular polarization, where the spin 
dynamics acquires a different sign. 

We gratefully acknowledge the financial support of 
the Deutsche Forschungsgemeinschaft through grant No. 
AX17/9-1. Financial support was also received from the 
Universidad de Buenos Aires, project UBACyT 2011- 
2014 No. 20020100100741, and from CONICET, project 
PIP 11220110100091. 


[1] E. Hemsing, M. Dunning, C. Hast, T. Raubenheimer, 
and D. Xiang, Phys. Rev. Lett. 113, 134803 (2014). 

[2] P. R. Ribic, D. Gauthier, and G. De Ninno, Phys. Rev. 



5 


Lett. 112 , 203602 (2014). 

[3] G. Gariepy, J. Leach, K. T. Kim, T. J. Hammond, 
E. Fmmker, R. W. Boyd, and P. B. Gorknm, Phys. Rev. 
Lett. 113 , 153901 (2014). 

[4] M. Padgett, Proceedings of the Royal Society of London 
A: Mathematical, Physical and Engineering Sciences 470 
(2014). 

[5] M.-E. Conprie, Nature Physics 9, 531532 (2013). 

[6] S. Patchkovskii and M. Spanner, Nature Physics 8, 
707708 (2012). 

[7] J. Wang, J.-Y. Yang, 1. M. Fazal, N. Ahmed, Y. Yan, 
H. Huang, Y. Ren, Y. Yue, S. Dolinar, M. Tur, and 
A. E. Willner, Nature Photonics 6, 488496 (2012). 

[8] N. V. Bloch, K. Shemer, A. Shapira, R. Shiloh, 1. Juwiler, 
and A. Arie, Phys. Rev. Lett. 108 , 233902 (2012). 

[9] P. Genevet, J. Lin, M. A. Kats, and F. Gapasso, Nature 
Gommunication 3, 1278 (2012). 

[10] U. D. Jentschura and V. G. Serbo, Phys. Rev. Lett. 106 , 
013001 (2011). 

[11] M. Padgett and R. Bowman, Nature Photonics 5, 343348 

( 2011 ). 

[12] B. J. McMorran, A. Agrawal, 1. M. Anderson, A. A. Herz- 
ing, H. J. Lezec, J. J. McClelland, and J. Unguris, Sci¬ 
ence 331 , 192 (2011). 

[13] J. Verbeeck, H. Tian, and P. Schattschneider, Nature 
467 , 301304 (2010). 

[14] M. Uchida and A. Tonomura, Nature 464, 737 (2010). 

[15] G. Molina-Terriza, J. P. Torres, and L. Torner, Nature 
Physics 3 , 305 (2007). 

[16] A. T. O’Neil, 1. MacVicar, L. Allen, and M. J. Padgett, 
Phys. Rev. Lett. 88, 053601 (2002). 

[17] H. He, M. E. J. Friese, N. R. Heckenberg, and 
H. Rubinsztein-Dunlop, Phys. Rev. Lett. 75, 826 (1995). 

[18] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and 
J. P. Woerdman, Phys. Rev. A 45, 8185 (1992). 

[19] G. F. Quinteiro and P. 1. Tamborenea, EPL 85, 47001 
(2009). 

[20] G. F. Quinteiro and P. 1. Tamborenea, Phys. Rev. B 79, 
155450 (2009). 

[21] G. F. Quinteiro and J. Berakdar, Opt. Express 17 , 20465 
(2009). 

[22] G. F. Quinteiro and P. 1. Tamborenea, Phys. Rev. B 82, 
125207 (2010). 

[23] G. F. Quinteiro, P. 1. Tamborenea, and J. Berakdar, 
Opt. Express 19 , 26733 (2011). 


[24] J. Watzel, A. S. Moskalenko, and J. Berakdar, Opt. 
Express 20 , 27792 (2012). 

[25] B. Sbierski, G. F. Quinteiro, and P. 1. Tamborenea, J. 
Phys.: Condens. Matter 25 , 385301 (2013). 

[26] M. B. Farias, G. F. Quinteiro, and P. 1. Tamborenea, 
Eur. Phys. J. B 86, 432 (2013). 

[27] G. F. Quinteiro and T. Kuhn, Phys. Rev. B 90, 115401 
(2014). 

[28] Y. Ueno, Y. Toda, S. Adachi, R. Morita, and T. Tawara, 
Opt. Express 17 , 20567 (2009). 

[29] N. B. Clayburn, J. L. McCarter, J. M. Dreiling, 
M. Poelker, D. M. Ryan, and T. J. Gay, Phys. Rev. 
B 87 , 035204 (2013). 

[30] F. Meier and B. P. Zakharchenya, eds.. Optical Orien¬ 
tation, Modern Problems in Condensed Matter Sciences, 
Vol. 8 (North-Holland Amsterdam, 1984). 

[31] C. Flindt, A. S. Sprensen, and K. Flensberg, Phys. Rev. 
Lett. 97 , 240501 (2006). 

[32] S. D. Ganichev, E. L. Ivchenko, S. N. Danilov, J. Eroms, 
W. Wegscheider, D. Weiss, and W. Prettl, Phys. Rev. 
Lett. 86, 4358 (2001). 

[33] P. Szumniak, S. Bednarek, B. Partoens, and F. M. 
Peelers, Phys. Rev. Lett. 109 , 107201 (2012). 

[34] 1. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76 , 323 (2004). 

[35] K. Schmalbuch, S. Gobbels, P. Schafers, C. Rodenbiicher, 
P. Schlammes, T. Schapers, M. Lepsa, G. Giintherodt, 
and B. Beschoten, Phys. Rev. Lett. 105 , 246603 (2010). 

[36] F. G. G. Hernandez, G. M. Gusev, and A. K. Bakarov, 
Phys. Rev. B 90 , 041302 (2014). 

[37] S. Kuhlen, K. Schmalbuch, M. Hagedorn, P. Schlammes, 
M. Patt, M. Lepsa, G. Giintherodt, and B. Beschoten, 
Phys. Rev. Lett. 109 , 146603 (2012). 

[38] E. Poem, O. Kenneth, Y. Kodriano, Y. Benny, S. Khat- 
sevich, J. E. Avron, and D. Gershoni, Phys. Rev. Lett. 
107 , 087401 (2011). 

[39] A. Greilich, S. E. Economou, S. Spatzek, D. R. Yakovlev, 
D. Reuter, A. D. Wieck, T. L. Reinecke, and M. Bayer, 
Nature Physics 5, 262 (2009). 

[40] Y. A. Bychkov and E. 1. Rashba, Journal of Physics C: 
Solid State Physics 17 , 6039 (1984). 

[41] G. Dresselhaus, Phys. Rev. 100 , 580 (1955). 

[42] G. Bastard, Wave Mechanics Applied to Semiconductor 
Heterostructures, Monographies de physique (Les Edi¬ 
tions de Physique, JOUVE, France, 1990) pp. 35-54. 



